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Abstract 

o , 

* ^ The CERN Axion Solar Telescope (CAST) experiment searches for axions from the Sun converted into photons with energies 
^^up to around 10 keV via the inverse Primakoff effect in the high magnetic field of a superconducting Large Hadron Collider 
t (LHC) prototype magnet. A backside illuminated pn-CCD detector in conjunction with an X-ray mirror optics is one of the three 
Ondetectors used in CAST to register the expected photon signal. Since this signal is very rare and different background components 
(environmental gamma radiation, cosmic rays, intrinsic radioactive impurities in the set-up, . . . ) entangle it, a detailed study of 
I the detector background has been undertaken with the aim to understand and further reduce the background level of the detector. 
The analysis is based on measured data taken during the Phase I of CAST and on Monte Carlo simulations of different background 
\Q components. This study will show that the observed background level (at a rate of (8.00 =b 0.07) x 10 -5 counts cm -2 s _1 keV -1 
between 1 and 7 keV) seems to be dominated by the external gamma background due to usual activities at the experimental site, 
while radioactive impurities in the detector itself and cosmic neutrons could make just smaller contribution. 
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1. Introduction 

Axions are a direct consequence of the Peccei-Quinn 
mechanism [1] proposed to solve the so-called strong CP 
problem (CP violation in strong interactions does not seem 
to exist in nature, although the QCD Lagrangian contains 
CP-violating terms). These particles could couple to two 
photons, which allows the production of axions inside the 
hot plasma of stars via the Primakoff effect. The expected 
flux of axions from the Sun has a mean energy of ~4.2 keV 
and virtually vanishes above 10 keV. In the presence of a 
transverse magnetic field, solar axions could be converted 
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back to observable photons via the inverse Primakoff effect, 
allowing axion detection on Earth [2] . 

The CERN Axion Solar Telescope (CAST) is intended to 
search for solar axions based on this principle [3,4]. It uses 
a decommissioned 9.26 m long LHC superconducting mag- 
net providing a 9T magnetic field, where the probability 
of axion-to-photon conversion is proportional to the length 
and magnetic field strength squared. Each bore of the twin 
aperture magnet has a cross-sectional area of 14.5 cm 2 . 
The magnet is installed on a platform that permits hori- 
zontal movement from azimuthal angle 46° to 133° and a 
vertical movement of ±8°. As a consequence, we can fol- 
low the Sun three hours per day, 1.5 hours each during 
sunset and during sunrise. The rest of the day is dedi- 
cated to background measurements. Three X-ray detectors 
are installed at each end of the magnet to search for ex- 
cess X-rays coming from axion conversions inside the mag- 
net during alignment with the Sun. The detectors are a 
Time Projection Chamber (TPC) [5] , a Micromesh Gaseous 





Fig. 1. Left: The pn-CCD detector with its vacuum housing. The pn-CCD chip is the central black part surrounded by a gold plated cooling 
mask. Towards the top the connection to a cold finger of a Stirling cooler device is provided by copper leads. The connection to the detector 
electronics is provided by a flex-lead (towards the bottom). The vacuum vessel is made of Aluminum. Right: The same picture with the 
internal passive copper-lead shield being installed (view from the X-ray optics). The circular hole in the shield provides the aperture towards 
the X-ray optics. The external lead shield is not shown. 



Structure (MICROMEGAS) [6], and a Charge Coupled De- 
vice (CCD) in combination with a X-ray mirror telescope 
[7]. The CCD and MICROMEGAS detectors observe sun- 
rise axions, while the TPC detector, covering two magnet 
bores, looks for sunset axions. During 2003 and 2004 the 
experiment operated with vacuum inside the magnet pipes 
(CAST Phase I). Due to coherence effects this setup allows 
us to explore an axion mass range up to 0.02 eV. No signal 
above background was observed in 2003, implying an upper 
limit to the axion-photon coupling of g a77 < 1.16x10" 10 
GeV -1 [8]. In order to extend the CAST sensitivity to 
higher axion masses, the CAST experimental setup has 
been transformed to be able to fill the axion conversion vol- 
ume with a buffer gas [4] , probing the axion mass range up 
to 0.8 eV (CAST Phase II). CAST started taking data in 
this configuration in November 2005. 

The CAST experiment is located at sea level at LHC 
point 8 of CERN and the investigated axion signal is ex- 
pected to be very rare. Therefore, different background 
components will be merged with the hypothetical signal in 
the detector data (environmental gamma radiation, cosmic 
rays, intrinsic radioactive impurities in the detector set-up, 
. . . ) . The goal of this work has been to study the back- 
ground of the pn-CCD detector of CAST with the aim to 
understand its origin and to further reduce the background 
level to increase the detector sensitivity. Data acquired dur- 
ing the 2003 and 2004 periods have been analyzed and a 
series of Monte Carlo simulations have been performed us- 
ing mainly the GEANT4 package [9] to achieve this goal. 
In particular, the response of this detector to photons and 
neutrons has been studied and a quantitative estimate of 
the contribution of different sources of background (like ex- 
ternal gamma and intrinsic radio-impurities in the detector 
materials) to the overall background of the CCD detector 
has been attempted. 



The structure of this paper is the following: details of 
the CCD detection system of CAST and its performance 
are presented in Sec. 2 and the main sources of background 
expected in CAST are described in Sec. 3, followed by a 
summary of the CCD background properties in Sec. 4. A 
description of the simulations and results is shown in Sec. 5. 



2. The CCD detector and the X-ray telescope 

The X-ray telescope of CAST is the most sensitive of the 
three different detector systems currently in operation, [4]. 
It consists of a Wolter I type X-ray mirror optics [10,11,12] 
focusing a potential axion signal on a small area on a CCD 
detector which is located in the focal plane of the optics. 
The CCD is a fully depleted back side illuminated pn-CCD 
with a depletion depth of 280/im and a pixel size of 150 x 
150 /im 2 , optimized for the 0.2-10 keV energy range and 
currently in use on board of ESA's XMM-Newton X-ray 
observatory [13]. The very thin SiC>2 radiation window of 
this detector with a thickness of only 30 nm, results in a 
quantum efficiency of « 1 between 1 and lOkeV. Such a 
device is operated for the first time in a low background 
experiment. For an in-depth description of the design of 
the X-ray telescope and its performance we refer to [7]. 
The long exposure time in CAST (> 3 years of operation 
time) provides a unique opportunity to study the long term 
background, performance, and different shielding concepts 
for such a detector in detail. 

The combination of a detector with high quantum effi- 
ciency and a focusing device already improves the signal to 
noise ratio of the experiment by a factor of « 12 compared 
to a non-focusing system. Further enhancement of the sen- 
sitivity could be achieved by adding a passive shield con- 
sisting of a combination of internal (inside the vacuum ves- 



sel) and external lead-copper components. As raw material 
for the shield components we have chosen low activity oxy- 
gen free copper and low activity lead, almost free of 210 Pb. 
In Fig. 1 the CCD detector is shown as it is integrated in 
its vacuum aluminum vessel with (right) and without (left) 
internal shielding components. 



3. Background sources 

The CAST experiment is located at one of the buildings 
of the point 8 experimental area at CERN. Therefore, any 
component of the environmental background at sea level 
can in principle affect the measured background of the de- 
tectors of CAST. 

The most important contribution is expected to be 
caused by external gamma rays, produced mainly by pri- 
mordial radio-nuclides like 40 K and the radioactive natural 
chains from 238 U, 235 U and 232 Th in laboratory soil, build- 
ing materials and experimental set-up as well as by 222 Rn 
in air. Gamma rays of cosmic origin make a negligible con- 
tribution. Activity levels of the walls of the experimental 
hall where CAST operates were measured using a germa- 
nium gamma spectrometer system [14] and radon levels 
have been also monitored there during long periods. Re- 
sults from these measurements will be taken into account 
later on to estimate the contribution of gamma rays to the 
total measured CCD background counting rate. 

Intrinsic radioactive impurities (either primordial or cos- 
mogenically induced) in the detector components and ma- 
terials can also make a relevant contribution in experiments 
looking for rare event signals because of their alpha, beta 
and gamma emissions. Since the CCD detector has an inter- 
nal shielding made of lead-copper within the vacuum cham- 
ber where it operates, impurities from the external compo- 
nents can in principle be disregarded since they should be 
greatly suppressed. Therefore, only impurities in the ma- 
terials composing the detector itself may be relevant and 
have been taken into account in this study. Activities from 
the CCD components were determined at the Canfranc Un- 
derground Laboratory using an ultra-low background ger- 
manium spectrometer and will be considered in evaluating 
their effect on the CCD measurements. 

Cosmic rays on the Earth's surface are dominated by 
muons and neutrons. While muon interactions (as those 
of other charged particles) can be rejected with « 100% 
efficiency thanks to their long ionizing tracks and the de- 
posited energy, signals from neutrons can contribute to the 
detector background. At sea level, most of the neutrons 
come from cosmic rays. The neutrons induced by cosmic 
muons or generated by fission or (alfa,n) reactions [15] are 
in principle marginal. However, since production of neu- 
trons by muons is strongly enhanced in high Z materials, 
muon-induced neutrons in the lead shielding have also been 
analyzed. 

In order to maximize the sensitivity of the CCD detector 
and to minimize the contribution of all background compo- 
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Fig. 2. A typical background spectrum measured with the CCD 
detector of CAST. Different fluorescent emission lines of material 
close to the CCD chip are apparent, e.g., the emission lines of copper, 
gold, and lead. 

nents, passive shields were installed and off-line rejection 
methods were applied. We want to point out that the signal- 
to-background ratio of the X-ray telescope is significantly 
better compared to the other two detectors in use in CAST 
due to the X-ray focusing system described in the previous 
section. Thus a sophisticated shielding concept including, 
e.g., a neutron shield seemed not to be of importance. 



4. Measured detector background 

During the last three years of detector operation in 
CAST, we have acquired more than 4298 h of background 
data with the X-ray telescope which have been partially 
analyzed. The 2003 background energy spectrum of the 
CCD detector, shown above in Figure 2, exhibits a count- 
ing rate of (11.5 ± 0.2) x 10 -5 counts cm -2 s _1 keV -1 . 
By adding internal and external shielding this level 
was reduced to a raw counting rate of (8.00 ± 0.07) x 
10 -5 counts cm -2 s _1 keV -1 . The spectrum shows several 
characteristic emission lines of the materials close to the 
CCD chip, e.g., the copper, gold, and lead lines, on top of 
an almost flat continuum spectrum (Compton scattered 
photons and secondary electrons [7]). In order to achieve 
a maximum sensitivity, the CCD detector is optimized for 
the energy range from 1 to 7keV which is the range of 
interest for the axion search. The major goal of this work 
is to explain the observed background in terms of different 
contributions from environmental background, neutron or 
muon induced background, and the contribution of natural 
radio-activity of the detector components themselves. 

Electronic noise of the detector and the so-called "par- 
tial" events have to be considered as intrinsic sources of 
background [16]. Since the electronic noise level of such 
CCDs is far below <C 100 eV, it is even for long exposure 
times not important for our energy region of interest with 
a minimum energy of E m - in > 0.7 keV. Calibrating spec- 
trum tell us that partial events contribute to it as a flat 



continuum shown in Fig. 5. These events arise from pho- 
tons absorbed in the detector near the Si-Si02 interface of 
the radiation entrance window (see reference [16]). For a 
mono-energetic line such events contribute to the observed 
spectrum with a constant fraction which is generally ex- 
pressed in terms of a peak-to- valley ratio (maximum of the 
photo-peak relative to the level of the flat-shelf). For the 
CAST pn-CCD detector this ratio is of the order of 2000. 
If we consider the Cu-Ka line of the background spectrum 
in Fig. 2 we would expect a contribution of at most « 
10 -6 counts s _1 keV -1 to the continuum spectrum between 
1-7 keV due to this effect. 

5. Background simulations 

In order to estimate the contribution of external back- 
ground and natural radioactivity to the overall background, 
a simulation tool for the CCD detector has been devel- 
oped using the GEANT4 package. As a first step we cross- 
checked the reliability of the code by comparing calibration 
data like, e.g., measured quantum efficiency and calibra- 
tion spectra with results we obtained from our simulations. 
The results will be presented in the following sections. Sim- 
ulations to understand the effect of the external environ- 
mental gamma background and its attenuation in different 
shielding configurations are presented and quantitative re- 
sults of the estimate of the contribution of this background 
to the CCD counting rates are shown. Results regarding 
neutrons, either from cosmic rays or muon-induced in lead, 
are also collected. Finally, an estimate of the contribution 
to the counting rate of the CCD detector of the radioactive 
impurities of the detector materials is made. 

5.1. The Code 

In the first simulations for neutrons and external gamma 
backgrounds a simplified description of the detector was im- 
plemented just considering the Si chip, the copper cooling 
mask and the printed circuit board. A much more detailed 
geometry for the detector was defined to carry out simula- 
tions of the radioactive impurities in the detector compo- 
nents; shapes and sizes of Si chip (including active area), 
ceramics, zero- force sockets, front and rear cooling mask 
(including gold lining), and printed circuit board have been 
reproduced as accurately as possible, keeping all the rela- 
tions between the sizes. Figure 3 shows views of the sim- 
ulated detector components. For both, the simplified and 
detailed implementations of the detector, it is placed in- 
side the aluminum housing. Copper-lead shielding as well 
as the tube connecting to the telescope have been consid- 
ered. Figure 4 shows the complete simulated set-up. 

Usual electromagnetic processes for photons and charged 
particles have been taken into account using models spe- 
cially developed for low energy (valid above 250 eV), which 
is important to reproduce fluorescence and emission of 
Auger electrons. Elastic and inelastic scattering, capture 
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Fig. 5. Measured (solid line) and simulated (dashed line) calibration 
spectra using a 55 Fe source. The escape peaks corresponding to the 
Cu-Koi and Cu-K/3 photo peaks are well reproduced in energy and 
integral intensity. The continuum apparent up to channel 125 is not 
reproduced by the simulation (see text). 

and fission are considered for neutrons. Radioactive de- 
cays are simulated also by GEANT4. Different simulation 
cut^f] have been defined for different kind of particles and 
regions, to save computing time while having accurate 
results. For each simulation, the spectrum of the energy 
absorbed in the CCD detector (considering the real inte- 
gration time of 72 ms) is recorded. 

5.2. Verification Against Calibration Measurements 

An energy calibration with a 55 Fe source was first sim- 
ulated and compared to the corresponding experimental 
spectrum. Figure 5 shows this comparison. The area of the 
peak at ^4.2 keV due to the escape of Si X-rays produced 
following a photoelectric interaction of 5.9 keV photons is 
quite well reproduced (differences are of 1.5%). The escape 
peak corresponding to 6.5 keV photons is probably embed- 
ded in the background and the other peaks. The continuum 
in the real calibration is much higher than in the simula- 
tion. This discrepancy is related to the partial charge collec- 
tion of some events, which affects the response of the CCD 
detector to low energy photons, as previously discussed in 
Sec. 4 (details about this can be found in Ref. [16]). This ef- 
fect accounts quantitatively for the observed ratio between 
the continuum level and the area of the biggest peak in cal- 
ibration measurements, but its inclusion in the simulation 
was disregarded since its influence on background measure- 
ments can be neglected. Other effects which could produce 
events in the continuum during calibrations (environmen- 
tal background, Compton continuum associated to peaks, 
energy losses of photons within the source collimator or in- 
ternal bremsstrahlung from the 55 Fe source) have been also 
studied, but their contribution would be much smaller. In 



2 Particles having an energy (or a path) lower than the fixed cut 
are no longer simulated and therefore they are directly absorbed. 



Fig. 3. View of all simulated CCD detector components together (left)and independent views of each of them (right): electronic board, 
ceramics, Si chip and front and back cooling mask. 




Fig. 4. Left: view of all simulated components, including shields, aluminum vessel and connection to telescope. Right: view of simulated 
components inside the aluminum vessel, including copper and lead shields and detector 
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Fig. 6. Simulated detection efficiency of the CCD to photons in the 
energy range from 0.2-100 keV (dashed line and triangles). In com- 
parison the measured quantum efficiency in the 1-15 keV is shown 
as well (solid line and diamonds). 

the measured spectrum a line from Al-K X-rays at ~ 1.5 keV 
is visible, due to the capsule where the calibration source 
is placed. That capsule with aluminum filter inside was not 
considered in the simulation. 

In addition to 55 Fe measurements, measurements to de- 
termine detection efficiency as a function of energy (see 
Ref. [13]) were also simulated and results are presented in 
Fig. 6. Independent simulations have been carried out with 
slightly different geometry details. Discrepancies between 
real data and simulations are of the order of ~4% below 
10 keV. Although no comparison with real data is possi- 
ble, high energy photons were also simulated obtaining for 
instance an efficiency as low as ~0.4% for 100 keV photons. 

5.3. External Gamma Background 

In order to evaluate the contribution of the environmen- 
tal gamma background to the CCD counting rate and the 
effect of the different shielding conditions, first the response 
of the detector has been studied by means of photon simu- 
lation. 

5.3.1. Response 

The response of the CCD detector has been studied by 
simulating photons with discrete energies from 30 up to 
2000 keV from the surface of a sphere containing all the 
components included in the simulation. Four different ge- 
ometries have been considered corresponding to the succes- 
sive shieldings added to the CCD set-up: no shield, inter- 
nal copper shield, internal copper shield and external lead 
shields and finally internal lead and copper shields and ex- 
ternal lead shield. Figure 7 presents the ratio between the 
counts in the region up to 7 keV (to exclude X-rays from 
copper) and the total simulated events for each one of the 
considered initial energies of the incident photons. Results 
for each one of the four different shielding conditions are 
plotted. 
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Fig. 7. Simulated response of the CCD detector to the external 
gamma background: events in the region up to 7 keV per sim- 
ulated photon as a function of the energy of the incident pho- 
tons. Four shielding conditions have been considered: with no shield 
(squares), with the internal copper shield (circles), with the inter- 
nal copper+external lead shields (triangles) and with internal cop- 
per+external lead+internal lead shields (crosses). 

Gamma rays with energies around 100 keV contribute 
most to the CCD background. However, the effect of this 
external gamma background above 50 keV, according to 
Fig. 7, does not seem to be very dependent on energy. This 
fact may be due to the balance of two trends: low energy 
photons are less penetrating but their contribution is higher 
while high energy photons can pass more easily through the 
housing and the shielding but their contribution is lower. 

The average reduction achieved in the counting rate be- 
low 7 keV due to environmental photons when using the 
complete shielding in the simulations seems to be slightly 
underestimated compared to the reduction actually ob- 
served in the CAST CCD detector (around a factor 3). Sim- 
plifications in the implementation of the components of the 
set-up, uncertainties in the response of the detector to high 
energy photons and/or non-isotropic background sources 
as considered in the simulation could be responsible for this 
discrepancy. 

5.3.2. Contribution to CCD counting rate 

Using this simulated response to external gamma includ- 
ing the effect of the complete shielding (see Fig. 7), an es- 
timate of the contribution to the counting rate of the CCD 
detector from the external gamma background has been at- 
tempted. Results of the radon level measurements made by 
the CAST Collaboration at the experimental site and the 
activity of the walls of the experimental hall measured us- 
ing a Ge spectrometer [14] have been used as input in this 
estimate. 

- Radon contribution: a mean radon concentration of ~10 
Bq m -3 has been assumed, according to available mea- 
surements. The rate of photons produced by 222 Rn de- 
scendants from the air in successive spherical layers of 
radius R has been evaluated up to R=10 m. From Fig. 7, 



the number of counts up to 7 keV produced by a photon 
of fixed energy coming from a spherical surface is known; 
this information has been used to derive the contribution 
to counting rates of radon emissions of different energies 
from each of the spherical layers described, just properly 
rescaling the sphere used in the simulations to the corre- 
sponding air spheres. Total contribution from radon up 
to 7 keV turns out to be ~ 10 -6 counts cm -2 s _1 keV -1 , 
almost negligible compared to the measured background 
rate in the CCD detector. 
- Walls activity contribution: specific activities for the 
North, South and East walls of the SR8 hall at CERN 
where the CAST experiment is running are known for 
40 K and for 232 Th, 235 U and 238 U chains (in some cases, 
since equilibrium is broken for 238 U chain two values are 
available up to and from 226 Ra) [14]. The lower part of 
the East and South walls (around 1/4 of the total height) 
is made of composite wall while the upper part is made 
of concrete, having different specific activities accord- 
ing to measurements. The North wall contains partially 
metal. The rate of photons coming from the surface walls 
from each of the measured radioactive impurities has 
been evaluated using these specific activities, the wall 
dimensions and a mean density of 2.3 g cm -3 . Uniform 
distribution of these impurities in the walls has been 
assumed and attenuation of photons through the walls 
has been taken into account. The contribution to the 
CCD counting rate of photons from each wall has been 
evaluated using again information from Fig. 7, rescaling 
the sphere used in the simulation to spheres having a 
radius given by the mean distance from the wall to the 
CCD detector. 

Table 1 summarizes the counting rates obtained as- 
suming the CCD detector at two different positions 
(height of 2.25 m, distance to South wall of 1.5 m and 
distance to East wall of 10 m for the first case and height 
of 3.70 m, distance to South wall of 0.5 m and distance 
to East wall of 9 m for the second one). Errors quoted 
in Table 1 (around 6%) come just from the error in the 
specific activities measurements. Considering the CCD 
detector at different positions along the tracking path 
when following the Sun does not change very signifi- 
cantly the contribution to the counting rates. A counting 
rate of around 3-4x 10 -5 counts cm -2 s _1 keV -1 could 
be attributed to the measured activities of the experi- 
mental site walls. This could therefore be the bulk of the 
measured background rate. 



5.4. Neutrons 

In a similar way as for photons, an evaluation of the con- 
tribution of neutrons from different sources to the CCD 
counting rate has been attempted, after studying the re- 
sponse of the detector to neutrons using simulations. 



Table 1 

Counting rates in units of 10 -7 counts cm -2 s _1 keV -1 expected 
from natural radioactivity of the walls of the experiment area for 
two different positions. 

Wall Height 2.25 m Height 3.70 m 

North 22.27 ±0.48 20.72 ± 0.45 

East (lower part) 3.92 ± 0.16 4.11 ± 0.17 

East (upper part) 45.8 ± 4.0 52.9 ± 4.6 

South (lower part) 21.96 ± 0.92 23.15 ± 0.97 

South (upper part) 187 ± 15 317 ± 26 

Total 281 ± 16 418 ± 27 



5.4.1. Response 

As for photons, the fraction of events detected in the 
CCD detector in the region of interest has been evaluated 
considering neutrons of different initial energies (from ther- 
mal up to 10 MeV). The different capacity of nuclei and 
electrons to ionize in the same medium compels us to take 
into account a relative efficiency factoiPH to compare en- 
ergy spectra produced by neutrons and by photons. An an- 
alytical expression based on the Lindhard theory and taken 
from Ref. [17] has been used for including this factor, which 
turns out to be ~0.28 for Si nuclear recoils of 10 keV and 
~0.36 for 100 keV. Nuclear recoils of up to ~32 keV will be 
registered in the region up to 10 keV of visible energy. 

The main contribution of neutrons to the counting rates 
in the energy region of interest comes from the elastic scat- 
tering off of silicon nuclei. Simple kinematics says that for 
neutrons up to ~250 keV all the induced recoils deposit in 
the detector a (visible) energy below 10 keV; neutrons with 
higher energies can produce more energetic recoils. 

Plots in Fig. 8 present the ratio between the events 
recorded in the region up to 10 keV of visible energy and 
the total simulated events, as a function of the energy of 
incident neutrons (up to 1 keV on top and above 1 keV at 
the bottom). The maxima appearing in Fig. 8 (bottom) are 
due to resonances in the cross-sections for neutron elastic 
scattering on silicon. Since the cosmic-ray neutron energy 
spectrum is peaked at thermal energies (~0.025 eV) [18], 
the response of the CCD detector to thermal neutrons 
has been also evaluated. The contribution of thermal neu- 
trons to the background in the energy region of interest 
(see Fig. 8, top) is at least one order of magnitude lower 
than the contribution of neutrons above 1 keV. Events 
induced by these very low energy neutrons are due not to 
nuclear recoils following elastic scattering but to radiative 
capture around and/or inside the detector; consequently, 
any gamma shield will help to further reduce the contri- 
bution of thermal neutrons. Below 1 keV, the contribution 
to background increases when energy decreases because of 
the higher relevance of neutron capture. 



3 It is defined as the ratio between the amplitudes of pulses produced 
by nuclear and electronic recoils which deposit the same energy 
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Fig. 8. Simulated response of the CCD to neutrons in the low energy- 
region: events registered up to 10 keV (visible) per simulated neutron 
as a function of the energy of the incident neutrons. On top, for 
neutrons up to 1 keV, at the bottom for neutrons above 1 keV. 

Table 2 

Measured neutron fluxes in units of 10 -3 cm -2 s -1 at sea level for 
different energy ranges. 



Energy Range Neutron Flux Reference 



< leV 


~ 1 


[18] 


0.4-0.1 MeV 


2.9 


[15] 


0.1-1 MeV 


1.6 


[15] 


1-10 MeV 


1.7 


[15] 


> 10 MeV 


5.6 


[19] 



5.4.2. Contribution to CCD counting rate 

A rough estimate of the contribution of the environmen- 
tal neutron background to the CCD counting rate in the 
region of interest has been attempted considering the neu- 
tron fluxes reported in Table 2 (corresponding to integrals 
in different regions of energy of the neutron flux at sea level 
according to Refs. [15,18,19]). Assuming for each one of the 
four energy intervals in Table 2 an average value for the 
fraction of events detected in the region of interest (from 
the response plotted in Fig. 8), the counting rate due to the 
external neutron background in the region up to 10 keV 
has been estimated to be 6 x 10 -6 counts cm -2 s _1 keV -1 . 
Since this number is one order of magnitude lower than 
the experimental background levels of the CCD detector, it 
seems that environmental neutrons are not very significant 
at the present level of sensitivity. 

At sea level, neutron production by capture of negative 
muons is strongly enhanced especially in high Z materials. 



Therefore, an estimate has been made also for the muon- 
induced neutrons in the lead shielding of the CCD detector. 
Induced neutrons come mainly from evaporation and there- 
fore their energy spectrum is peaked around 1 MeV and 
very reduced above 5 MeV. Using a measured total muon 
flux of - 50.3 x 10" 3 cm -2 s" 1 , FLUKA simulations [20] [21] 
give a yield for neutrons of around ~ 1 x 10 _3 cm _2 s _1 
which is more than one order of magnitude lower than the 
flux due to environmental neutrons. Therefore, the instal- 
lation of the present lead shield is not a problem from the 
point of view of muon-induced neutrons. 

It is worth noting that the values quoted for the esti- 
mated contribution to the CCD counting rate of the exter- 
nal gamma and neutron backgrounds must be considered 
as just an indication of the order of magnitude of these 
contributions. Both particle fluxes and simulated detector 
responses (the key ingredients of these estimates) have im- 
portant uncertainties which in many cases are difficult to 
be quantified. Some facts that contribute to these uncer- 
tainties are that high energy photon and neutron simu- 
lations have not been checked against experimental data, 
that radon levels are known to have important fluctuations 
in time or that the neutron flux at sea level depends on dif- 
ferent factors not taken into account. Nevertheless, results 
have been useful to identify the main background sources 
of the CCD detector. 



5.5. Intrinsic radioactive impurities 

An evaluation of the contribution of the internal radioac- 
tive impurities in the components of the CCD to the count- 
ing rate measured in this detector has been made based on 
some measurements of radiopurity levels of several compo- 
nents and simulation of the effect of these activities. These 
estimates can help in the development and construction of 
a new CCD detector optimized from the radioactive point 
of view. 

The levels of radioactive impurities in the main compo- 
nents of a CCD detector equivalent to that actually used in 
the CAST experiment (circuit board, sockets, ceramics, Si 
detector and cooling mask) were measured in the Canfranc 
Underground Laboratory in Spain using a special ultra- 
low background germanium detector and can be found in 
a database of radiopurity of material£l]inside the ILIAS 
program (Integrated Large Infraestructures for Atroparti- 
cle Science). Activities come mainly from the radioactive 
chains 235 U, 238 U, 232 Th and the isotope 40 K. For the Si 
chip and the cooling mask pieces only upper limits could 
be derived for the impurities of the radioactive chains 235 U, 
238 U and 232 Th. For the chip, also for 40 K only a limit was 
deduced. 238 U chain is broken at the level of 226 Ra for the 
sockets and ceramic pieces. 



4 Available at |http : / /radiopurity . in2p 3 . f rj] 
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Fig. 9. Simulated energy spectrum registered in the CCD for 232 Th 
impurities in the ceramics. 

5.5.1. Simulations 

For each component of the CCD detector, complete 
decays of the radioactive chains 235 U, 238 U and 232 Th, 
assumed in secular equilibrium, and the isotope 40 K 
have been simulated considering impurities uniformly dis- 
tributed in the materials. As a result, a collection of energy 
spectra corresponding to each impurity in each component 
has been obtained. Just as an example, Fig. 9 shows the 
simulated spectrum registered in the CCD detector for 
232 Th impurities in the ceramic pieces (it must be noted 
that no normalization has been made in this plot). 

Peak structures observed in the real CCD spectrum (see 
Fig. 2) at energies slightly higher than those of the region 
of interest are partly reproduced by these simulations: 

- In most of the simulated spectra for the different compo- 
nents, the X-ray K a line from copper at 8.0 keV is clearly 
imposed on a flat background. It is due to photoelectric 
interactions of more energetic photons in copper pieces. 
A much smaller line at 8.9 keV is visible too due to much 
less intense copper Kp X-rays. 

- X-rays from gold (L a at 9.7 keV and Lp at 11.5 keV) 
are also present since the gold lining of the cooling mask 
pieces where they are generated has been properly in- 
cluded in the simulations. 

- Other low energy gamma lines and especially X-rays be- 
tween 10 and 15 keV emitted by daughter nuclei in the 
radioactive chains have been identified in the simulated 
spectra for the components closest to the Si chip. But 
the X-rays observed in the measured CCD spectrum and 
presumably attributed to ordinary solders in the elec- 
tronic board containing 210 Pb are not reproduced since 
these solders have not been taken into account in the 
simulations. 

- In the simulated spectrum for 40 K in the chip, a peak at 
~3 keV is seen due to the X-rays following the electronic 
capture in this isotope. This peak is not observed in the 
real spectrum because it is probably hidden by the con- 
tinuum background, which points to a completely negli- 
gible 40 K impurity in the Si chip. 



Table 3 

Total contribution to CCD counting rate between 1 and 7 keV from 
natural radioactivity of the individual detector components in units 
of counts cm -2 s _1 keV _1 . We considered the 238 U, 235 U, 232 Th, 
and 40 K decay chains. 

Detector Component Differential Flux 



CCD Board 
CCD Chip 
Ceramics 
Sockets 
Front Cooling Mask 
Rear Cooling Mask 



(5.83 ±0.41) x 1CT 7 

< 2.20 x 10- 5 
(1.17 ±0.14) x 10- 6 
(1.34 ±0.15) x KT 6 

< 3.60 x 10- 7 

< 2.61 x 10- 7 



5.5.2. Results 

Table 3 shows the total contributions to the CCD count- 
ing rate of each simulated component between 1 and 7 keV. 
Contribution from the Si chip and the cooling mask pieces 
must be considered just an upper limit, since only bounds 
for activity were set for these components, as stated before. 
Although the screening of radiopurity was made only for 
the front part of the cooling mask, equal activity has been 
considered for the rear part assuming that the same kind 
of copper was used for both pieces. Contribution from 238 U 
in the sockets and ceramic pieces is overestimated since the 
assumed secular equilibrium in the chain is not true, hav- 
ing from 226 Ra onwards a lower activity. The errors quoted 
in Table 3 (between 5 and 10% approximately) include the 
statistical error associated with the simulation and also the 
error of the measured activity of the pieces. In general, er- 
rors from the activity dominate, although statistical errors 
are not always negligible because of the large execution 
time of simulations. 

These estimates indicate that the components whose 
radioactive impurities have been actually measured (CCD 
board, ceramics and sockets) produce (3.05 ± 0.20) x 



10" 



keV in the region up to 7 keV. Con- 



sidering the upper limits of radio-impurities for all the other 
components, the total contribution to background would 
be 2.6 x 10 -5 counts cm -2 s _1 keV -1 . Compared to the 
measured background level in the CCD detector between 
1 and 7keV, (8.00 ± 0.07) x 10~ 5 counts cm" 2 s" 1 keV~\ 
impurities from CCD board, ceramics, and sockets account 
for of the observed counting rate; including all the 

other simulated impurities, up to ~33% of the counting 
rate would be justified at most. 

6. Conclusions 

The CCD detector of the CAST experiment looks for 
the very rare signal of solar axions, consisting of photons 
peaking around 4 keV. Different background components 
entangle the expected signal. A Monte- Carlo simulation 
tool for this detector has been developed and used for the 
most relevant background sources with the aim to help in 
the understanding of the origin of observed events and con- 
structing a plausible background model for the measured 



counting rate (8.00 ± 0.07) x 10" 5 counts cm -2 s _1 keV -1 
between 1 and 7keV. 

Using the simulated response of the CCD detector to 
the external gamma background (including the effect of 
shieldings), an estimate of the contribution of this back- 
ground component has been attempted, finding that mea- 
sured radon levels in the air of the CAST site could pro- 
duce on average ~1% of the registered counting rate, while 
however just the measured activities from 40 K and 232 Th 
and 238 U chains in the walls of the CAST hall could justify 
more than 50% of this counting rate. 

The response of the CCD to neutrons of different energies 
has also been simulated. Based on this response and typi- 
cal fluxes of neutrons at sea level, a rough estimate of the 
contribution to the CCD counting rate of environmental 
neutrons has been made. At the present level of sensitivity, 
these neutrons do not seem to be a very relevant source of 
the CCD background, producing just a few per cent of the 
observed counting rate. Contribution from muon-induced 
neutrons in the present lead shieldings has been checked to 
be negligible. 

Finally, the contribution to the CCD counting rate of the 
internal radioactive impurities of the main detector compo- 
nents has been simulated using the activities measured at 
the Canfranc Underground Laboratory with an ultra-low 
background germanium detector. This contribution could 
justify at most 33% of the measured counting rates. 

Taking into account all these results, a quite complete 
model for the background measured by the CCD detector 
has been obtained. Other possible relevant sources of back- 
ground, not evaluated up to now, are thought to be the 
radioactive impurities from the soil and from massive com- 
ponents of the experimental set-up. 
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